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Abstract 

We present a theoretical study of the quantum magnetocapacitance of an ultrathin topological 
insulator film in an external magnetic field. The study is undertaken to investigate the interplay 
of the Zeeman interaction with the hybridization between the upper and lower surfaces of the thin 
film. Determining the density of states, we find that the electron-hole symmetry is broken when 
the Zeeman and hybridization energies are varied relative to each other. This leads to a change 
in the character of the magnetocapacitance at the charge neutrality point. We further show that 
in the presence of both Zeeman interaction and hybridization the magnetocapacitance exhibits 
beating at low and splitting of the Shubnikov de Haas oscillations at high perpendicular magnetic 
field. In addition, we address the crossover from perpendicular to parallel magnetic field and find 
consistency with recent experimental data. 
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Topological insulators (TIs) form a new class of compounds with strong spin-orbit in- 
teraction [l|-l8|. These materials possess an energy gap in the bulk (like insulators) and 
conducting electronic states at the surface (like metals). The nature of the surface states 
depends on the dimensions of the TI j9-16|: Two-dimensional surface states have been 



revealed in transport experiments for three-dimensional TIs. The charge carriers in such 
surface states form a two-dimensional gas of Dirac fermions with a single Dirac cone energy 
spectrum. Dirac fermions have been observed in angle resolved photo-emission experiments 
and quantization of their Landau levels (LLs) has been confirmed by scanning tunneling 



spectroscopy 



17|-2ll- 

TI thin films provide a new degree of freedom which is the thickness 



of the film. Effects of the crossover from 3D to 2D topological surface states as well as 
hybridization of the states at the two surfaces of the thin film have been studied for varying 
: ilm thicknesses and found to become important when the thickness is less than about 5 nm 
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26| . The Dirac fermions in TIs share several features with those of graphene 27|, |28 |. 



However, there are important differences between graphene and the surface states of a TI. 
Unlike graphene, there is an odd number of not spin degenerate Dirac points at any surface. 
In addition, TIs exhibit a unique spin-momentum locking in the surface states. 

Observation of magnetic quantum oscillations, Shubnikov de Haas (SdH) and de Haas 



van Alphen 



29 



30l |. in the transport experiments has resulted in theoretical work on mag- 



netization [29| and Berry phase 30|, |3l| features of TI surface Dirac fermions. Furthermore, 



much attention is paid to magneticall y d oped TIs due to their exceptional properties and 



potential spintronics applications [lo|, |24]. In this context, capacitance measurements are 
among the most important tools for studying the electronic properties. They can be effec- 
tively used to probe the thermodynamic density of states (DOS). Although the focus in TI 
research has been on transport properties, insight into the fundamental electronic properties 
and device physics still calls for knowledge about the capacitance-voltage (C-V) characteris- 
tics of the system. Results have been reported for carbon nanotubes, graphene nanoribbons. 



and mono- and bilayer graphene systems 



37|. In view of this, attention is now being 



paid to electrostatic properties, such as the magnetocapacitance of TI thin film devices. 
Furthermore, for improve the performance of field effect transistors jssi, the potential of TI 
thin films as channel materials is creating much excitement. This is due to the excellent 
intrinsic transport features as well as the possibility of patterning device structures within 
top-down lithographical techniques. The present work aims at determining the combined 



effects of the Zeeman interaction and hybridization on the magnetocapacitance of a TI thin 
film device. 



I. RESULTS 

We consider Dirac fermions on the surface of a 3D TI thin film which we take to be in 
the xy-plane. There is an external magnetic field perpendicular to the thin film and we 
take into account hybridization between the upper and lower surfaces. The two-dimensional 



Hamiltonian for Dirac fermions in a magnetic field with hybridization is given by 
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H = via^TTy - T^cTyTz^) + (r^A^ + Ah)(T^. (1) 

Here, and ay are Pauli matrices that operate on a combined space of spin and surface 
degrees of freedom, = +/— denotes symmetric/antisymmetric combinations of the upper 
and lower surfaces of the TI, v denotes the Fermi velocity of the Dirac fermions, tt = p+eA/c 
is the two-dimensional canonical momentum with vector potential A, and c is the speed of 
light. Moreover, we have the Zeeman energy = ^gjisB, the effective Lande factor 
g, the Bohr magneton /i^, and the hybridization matrix element A/^, which reflects the 
hybridization between the upper and lower surfaces of the TI. We employ the Landau gauge 
and express the vector potential as A = (0, Bx^ 0). The LL energies then are given by 

El^ = -{A, + T,/\h) (2) 



with A = +/— for electron/hole bands. Moreover, u = v\/eB/h is the cyclotron frequency 
of the Dirac fermions and n is an integer. The main point to note about the energy spec- 
trum, which will be crucial to the following discussion, is: The splitting of the energy levels 
depends on both the Zeeman interaction and the hybridization between the upper and lower 
surfaces. If either of these two mechanisms is absent, no splitting will occur. Further, the 
energy spectrum is electron-hole symmetric in the absence of the Zeeman interaction, with 
a hybridization gap. We note that angular resolved photoelectron spectroscopy finds no 
electron-hole symmetry even in the absence of both Zeeman interaction and hybridization, 
reflecting a quadratic term in the Hamiltonian 39|, |40|. However, in our surface Hamiltonian 



in Eq. (2) no quadratic term is present as it can be neglected when the system is doped such 



that the Dirac point comes close to the charge neutrahty point, which is the focus of our 
work. In the presence of both Zeeman interaction and hybridization, the n = LL exphcitly 
shows a quantum phase transition. This will be important in the discussion of the DOS and 
the capacitance results. 

We consider a gated TI device in which the capacitor is formed between the gate and the 
TI thin film. The magnetocapacitance of the device represents the charge response in the 
channel as the channel potential is varied. In conventional devices, the magnetocapacitance 
Cq is large and can be ignored. On the contrary, in low dimensional devices, such as 2D elec- 
tron gases, graphene, and future TI devices, it is the dominant capacitive contribution and 
therefore an important quantity in the design of nanoscale devices. The central expression 



for the quantum capacitance is 
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where ng is the carrier concentration and ep the Fermi energy. The temperature dependent 
DOS at a finite magnetic field, Dx{B), is obtained from the relation 



dfie-ep) . 
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with the Fermi Dirac distribution /. In the limit of zero temperature we have Dt^B) 
D^sp) and 
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Assuming a Gaussian broadening of the LLs, the DOS per unit area is given by 
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where F is the width of the Gaussian distribution (zero shift). For n = we have 
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where Eq = — (A^ + T^Ah). Combining Eqs. (8) and (9), we obtain 
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At the charge neutrahty point {ep = 0) the result is 
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which can be written as 
D{ep = 0) 



exp 
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where x = ft^i^^/^T"^. Equation (10) reduces to the result reported in Ref. [36| in the absence 
of Zeeman interaction and hybridization. 



The magnetocapacitance Cq is plotted in Figs. 1 and 2 as a function of t 
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le Fermi energy 
lll-ksl: g = 60, 



(i.e., of the gate voltage). The following parameters are employed 
i? = 3 T, A^ = 5 meV, and A^ = 3 meV. To obtain analytical results, we choose a constant 
level width of F = 0.3 meV. We are interested in changes of the character at = on 
variation of the Zeeman interaction and hybridization relative to each other. It must be 
noted that the n = LL here plays the most important role. Figure 1 shows that Cq is 
zero at ep = 0. This occurs because the n = LL splits only when the hybridization does 
not vanish. The n = LL splits into one electron and one hole level, which reflects a metal 
to insulator transition caused by the hybridization (for A^ = 0). The states are electron- 
hole symmetric at this stage. Cq can be tuned from a minimum to a maximum when the 
Zeeman energy is increased by changing the external magnetic field. For A^ > Ah both 
n = sublevels are located in the hole region, which breaks the electron-hole symmetry, see 
Fig. 2, and represents the trivial to non-trivial topological insulator phase transition. To 
observe the splitting, the broadening of the LLs must be less than the hybridization energy. 
We note that not only the n = LL but all LLs are split into two sublevels. 

In Fig. 3 we address the SdH oscillations in Cq. The following parameters are used: 
F = 0.3 meV, T = K, ne = 4 x 10^^ m'^, v = 3 x 10^ ms-\ and Ah = 4 meV. For low 
and high magnetic field we observe a beating pattern and a splitting of the SdH oscillations, 
respectively. Both occurs due to interference of the SdH oscillations at two different frequen- 
cies and therefore is a consequence of the splitting of the LLs for finite Zeeman interaction 
and hybridization. The beating pattern vanishes once the Zeeman interaction dominates 
the hybridization. For the chosen parameters this occurs at a magnetic field of about 1.5 
T. Above this value we find a well resolved splitting of the SdH oscillations in Fig. 3. A 



similar beating pattern and splitting of SdH oscillations in the longitudinal resistivity is seen 



in Ref. 4J], where it occurs because of splitting of the LLs as a result of a broken inversion 
symmetry due to substrate effects. The sample investigated in this reference is thick and 
hybridization effects are ruled out. 

In order to understand the origin of the beating pattern and splitting of the SdH oscilla- 
tions quantatively, we employ the Poisson summation formula 



oo „ oo „ 

-F(0) + ^F{k) = / F{n)dn + 2^(-l)'^ / F{n) cos[27Tkn\d- 

k=l n k=l n 



tn (11) 
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to arrive at the following expression for the DOS at zero temperature which incorporates 
the contributions of the n > 1 LLs: 

2' 
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Here, Dq = lepl/'n'fi^v'^ is the zero magnetic field DOS. For F ^ hw, it is sufficient to retain 
only the first order term {k = 1), since the higher order terms are highly damped. From Eq. 
(12) we can confirm that the beating pattern is due to the interference of waves with two 
different frequencies. The oscillatory part of the k = 1 term in the sum for = +/— is 



cos 
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and can be expressed as 
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This represents a wave of higher frequency whose amplitude oscillates at a lower frequency, 
giving rise to the beating pattern, because A^A/i « ep- The amplitude of the SdH oscil- 
lations is modulated by cos [^rl^A^A/i] and nodes occur at = ±0.5, ±1.5,..., where 
the modulating cosine vanishes. We also note that the amplitude modulation given by the 
term cos [^^I^A^A^] occurs when both Zeeman interaction and hybridization are present in 
the system. Moreover, the threshold magnetic field where beating is seen depends on the 
strenghts of the Zeeman interaction and hybridization. For magnetic fields down to about 
1.5 T the beating persists. Above this value it is quenched to show splitting of the SdH 
oscillations. As mentioned before this finding is in agreement with the experimental work 



in Ref. 



44|. 



II. DISCUSSION 



Topological systems are typically doped away from the charge neutrality point. In addi- 



tion, the bulk usually dominates the surface contributions. However^ 



insulators make it possible to access the surface contributions only 
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lin films of topological 



4l|-|43j , typically when 



the thickness is below 5 nm. In this context, we compare our results with the experimental 



observations for the quantum magnetocapacitance [45|] . The latter authors study films that 
are thick enough that hybridization between the top and bottom surfaces can be ignored (10 
nm Bi2Se3 samples). To compare their results for a tilted magnetic field, we next consider 
this situation and explore the crossover from external perpendicular to parallel (in-plane) 
magnetic field on the quantum capacitance. For a tilt angle 0° < < 90° with the normal 
of the film, the magnetic field is given by (i?sin6', 0, Bcos6). To describe the effects of the 
tilting we start from the two-dimensional Hamiltonian 

H = VT^ia^TTy - CXyTTr,) (15) 

where A = (0, xB cos 6 — zB sin 6, 0) and z is the thickness of the film. The LL energies are 
obtained as 

En,\ = X\/2nfvjjf, (16) 

with ujt = v\/ eB cos6/h. Note the dependence of the LLs on the tilt angle 6 as compared 
to Eq. (2). 

Following the same procedure as before, we employ the energy spectrum to compute 
the quantum capacitance, which is plotted in Fig. 4 as a function of the Fermi energy 
(gate voltage) for different tilt angles. We focus on the broadening of the LLs and the 
corresponding effects on the SdH oscillations. We fix i? = 7 T. The SdH oscillations are 
suppressed as we increase the tilt angle, since we are decreasing the perpendicular component 
of the magnetic field. For 6 — t- 90°, when the magnetic field is almost completely aligned 
with the film, the SdH oscillations are washed out. The broadening of the LLs generally 
depends on the magnetic field strength, LL index, and scattering parameters. This requires 
a self consistent calculation, which usually is performed numerically. In order to carry out 
a tractable analytical calculation we choose a constant level width of F = 25 K. The results 



in Fig. 4 are consistent with the experimental Fig. 4 in Ref. 45 1. 



An experimental investigation of transport arising from the surface states in TIs is hin- 



dered by the bulk contributions due to naturally occurring defects and residual carrier doping 



46N48| ■ In this regard, the low mobility of the surface states is a serious problem. However, 



recent work on TIs has shown a way to overcome this problem by studying transport in 
ultra thin samples, where a 12 times higher surface mobility as compared to the bulk can 
be achieved |49i]. It has been demonstrated that the large surface to volume ratio of Bi2Se3 
thin films effectively sup presses bulk effects and enables probing of the electronic transport 
of the surface states [SOf] . There has been a rapid improvement in the sample quality coming 
along with a reduction of film thicknesses 5l|. Relevant to the present study, the effect of 



the Zeeman exchange field in TIs is significant 
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51| with a magnitude around 50 meV due 



to the large effective g-factor. Also, recent experiments have found a hybridization gap in 



the ultrathin regime (less than 6 nm) 22|, |26| as well as quantum oscillations [43]. Therefore, 



it can be expected that the effects discussed before can be realized in quantum capacitance 
experiments in ultrathin films. 



In conclusion, we have investigated the effects of Zeeman interaction and hybridization 
on the magnetocapacitance of a TI thin film in a perpendicular magnetic field. To this 
aim, we have obtained an analytic expression for the DOS incorporating both the Zeeman 
and hybridization contributions. It turns out that a combination of Zeeman interaction and 
hybridization breaks the electron-hole symmetry, resulting in significant effects on the mag- 
netocapacitance at Ep = 0. In addition, we have shown that in the presence of both Zeeman 
interaction and hybridization, the magnetocapacitance exhibits beating at low and level 
splitting at high perpendicular magnetic field. This behaviour is explained by interference 
of the SdH oscillations at the two frequencies induced by the splitting of the LLs. Let us 
comment on the relevance of our findings for experiments. The cyclotron energy for i? = 1 
T is /iw = 36 meV. Observation of the splitting of the LLs and the discussed consequences 
requires that the temperature is low and the system clean enough that disorder effects do 
not wash out the splitting. Temperature and disorder broadening thus must not reach the 
Zeeman or hybridization energies. For a tilted magnetic field with 6 — t- 90° we find that the 
SdH oscillations are washed out. 
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FIG. 1. Quantum capacitance as a function of the Fermi energy atT = 0K,i? = 3T, and 
Zeeman energy meV. The hybridization energy is meV (dotted lines) and 4 meV (solid lines), 
respectively. 
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FIG. 2. Quantum capacitance as a function of the Fermi energy at T = K, B = 3 T, Zeeman 
energy 5 meV, and hybridization energy 3 meV. 

FIG. 3. Quantum capacitance as a function of the magnetic field at T = K and hybridization 
energy 4 meV. 

FIG. 4. Quantum capacitance as a function of the Fermi energy for different tilt angles of the 
magnetic field: 0° (dashed lines), 40° (red solid lines), 70° (blue solid lines), and 80° (black solid 
lines). We use 5 = 7 T and F = 25 K. 
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